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ABSTRACT 


The structural differences of the ion and neutral composition in 
the thermospheric region are studied by solving a system of basic iono- 
spheric and atmospheric equations. The study shows that the compositional 
changes during a magnetic storm arise largely as a result of changes in 
the neutral composition at the turbopause. A decrease in [0]/[Np] in 
the lower atmosphere triggers a complex chain of events which Results in 
an increase of the neutral gas temperature, depletion of the 0 layer and 
enhancement of N0 f . The relative changes in these layers occasionally 
produce a sequence of electron density profiles giving rise to the so- 
called G condition. It is shown that, compared to the neutral atmosphere, 
the ionosphere is mucn more sensitive to the changes in [0]/[N 2 ] in the 
lower thermospheric region. Since the ionospheric parameters can be 
measured much more accurately than the atmospheric parameters, it is 
argued that they should form an integral part of the observational data 
required to construct the atmospheric models. 


1 . INTRODUCTION 


It has been recognized for some time that the ionospheric effects 
associated with geomagnetic disturbances are principally the result of 
changes in the neutral composition of the upper atmosphere (King, 196?, 

1966, 1967; Duncan, 1969). This assumption has generally been invoked to 
explain the changes in the F-region critical frequencies (f Fg) observed 
during the storm period. Since direct measurenents of the neutral com- 
position during geomagnetic disturbances are almost non existent, the 
validity of this concept could not be tested beyond its ability to explain 
the observed increase or decrease in f F 0 . In recent years electron den- 
sity and temperature measurements from satellites and incoherent radar 
backscatter during magnetic storms nave extended our horizon beyond the 
height of the F-region maximum (King et al., 1967; Bauer and Kri shnamurthy , 
1968; Norton, 1969; Reddy et al . , 1967, 1969; Evans, 1965; Rao, 1968). 
Although these observations are only few, they provide a much wider base 
for testing the concept of neutral composition changes during geomagnetic 
storms. In addition, a number of important papers have appeared in recent 
years describing the behavior of airglow in the various wave length bands 
during magnetic storms (Weil and Glaume, 1967; Truttse, 1968). These 
measurements provide direct evidence to the changes in the neutral composi- 
tion (Krassovsky, 1968) and thus offer an excellent opportunity for studying 
the behavior of the neutral atmosphere during storm periods. 

The physical mechanisms causing variations in the neutral atmosphere 
during storms are not fully understood. Cole (1962, 1965) has proposed 
two principal mechanisms by which energy is deposited in the upper atmosphere 
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during magnetic storms. These are: (1) the thermal conduction of the 
energy from the magnetosphere to the ionosphere and; (2) the Joule dissi- 
pation of the ionospheric currents. The two mechanisms are complementary 
and not mutually exclusive. They are responsible for a deposition of 
energy in different regions of the upper atmosphere. Since Joule heating 
occurs principally in the auroral regions, the transfer of energy to the 
lower latitudes probably occurs by gravity waves (Hines, 1965). Cummings 
and Dessler (1967) propose Joule dissipation of the asymmetric ring 
current as an alternative mechanism for the deposition of energy in the 
auroral and subauroral region of the lower atmosphere. Regardless of the 
merit of any particular mechanism, it is reasonable to believe that the 
energy deposited in the altitude range around 120 km should cause a change 
in the circulation pattern of the lower atmosphere resulting in the 
change of the turbopause level and the neutral compositon of that region. 
This is the basic assumption of this paper. Its justification lies in 
its ability to explain the storm related changes in the upper atmosphere. 

Any change in the composition of the neutral atmosphere not only 
affects the structural parameters of the region where the change takes 
place but also affects the other regions profoundly through a series of 
complex processes involving diffusion, transport and heat conduction. 

Thus, in order to understand the changes in the upper atmosphere, it is 
important to consider the behavior of the neutral and ionized constituents 
as one problem. This point was clearly established in a recent paper by 
Chandra and Herman (1969). Through simultaneous steady state solutions 
of the electron continuity equation and the heat conduction equations for 
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electrons, ions and the neutral gas, the authors showed how some of the 
well known features of the F-region storm can be explained by assuming 
a decrease in the ratio of atomic oxygen to molecular nitrogen, [0]/ [N^] * 
at the lower boundary. Such an assumption was consistent with the storm- 
time increase in the neutral density inferred from satellite drag observa- 
tions . 

The purpose of this paper is to present detailed time dependent 
calculations of the changes in composition of both the ionized and neutral 
constituents of the upper atmosphere resulting from the changes in the rela- 
tive values of [0], [Og] and [N^] at the lower boundary. The results are 
discussed in terms of the observed changes during geomagnetic storms. 

2. FUNDAMENTAL EQUATIONS 

The results of this paper are based on the simultaneous solution of 
a system of coupled, time dependent differential equations consisting of: 

1. Equations of continuity for 0 + , N0 + , and H + 

2. Heat conduction equations for 0 + , H + and the electron and neutral 
gas. 

3. Equations of motion of the various ionic and neutral constituents. 

These equations and their method of solution have been fully described 

by Stubbe (1970) and will not be repeated here. The neutral atmosphere is 
assumed to consist of 0^, N^, 0, H and He. Their height variaions are 
described by Chandra and Stubbe (1970) in terms of the dynamic diffusion 
model. The dynamic diffusion model was introduced in order to account for 
the observed phase difference between the neutral density as derived from 
satellite drag data and the neutral temperature as derived from incoherent 
backscatter. 
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A J the calculations are performed for a mid-latitude station (Geo- 
graphic latitude 42°; dip angle 72°;, roughly corresponding to Millstone, 
at equinox. The solar flux is assumed to be 30% higher than that given 
by Hinteregger et al . ( 19 65) » corresponding to low to medium solar activity. 
Following the convention of Chandra and Herman, we have characterized the 
assumed values of [0], [^L [N ? ] and at the lower boundary (120 km) 
by quiet and disturbed conditions. Their values at 120 km for 1800 hours 
local time are given in the following table. 



QUIET (Q) 

DISTURBED (0,) 

DISTURBED (D 2 ) 

C\J 

0 

1 .1 3xl0“/cm 3 

1 .13xl0“/cm 3 

1 .69x10* ‘/cm 3 

N 2 

6 . 00 x 1 0 1 ‘/cm 3 

6 . 00 x 10 * ‘/cm 3 

9.00X10 1 Vcm 3 

0 

6.76x10* °/cm 3 

3.80xl0‘°£m 3 

5.00x10* °/fcm 3 

t n 

355°K 

355°K 

355°K 


The time variations in density at 120 km are assumed to be of the form 
(1+0.3 coswt), and the temperature is assumed to vary as (1+0.1 cosait) where 
00 is the angular frequency of the earth's rotation and t the time measured 
from noon. Note that the difference between D-j and Q is only with respect 
to the values of [0]. In D-j , the [0] concentration is assumed to be about 
50% less than in Q. In D^, [0] is 25% less than in Q, while [ 03 ] and [N^] 
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are assumed to be increased by 50%. 

The composition changes from Q to and conditions assumed in 
this paper are arbitrary and are meant for illustration only. The obser- 
vational data on the neutral composition in the lower atmosphere is still 
scarce. The difference between the various measurements is so large that 
it is difficult to construct models for quiet conditions, much less for 
disturbed conditions. The Tasks of making quantitative estimates for com- 
position changes during the period of geomagnetic storms is even more for- 
midable. The altitude of the turbopause level and its variation during 
magnetic storms are the key parameters which must be calculated to compute 
the changes in the neutral composition at a specified reference level. 

This requires not only the knowledge of the source of energy and its mag- 
nitude but also the mechanism for its dissipation before the system of 
dynamical and energy equations can be solved for this purpose. Evidence 
has recently been presented (Testud and Vasseur 1970) that the energy 
generated in the auroral thermosphere during magnetic storms can be carried 
to the lower latitudes via gravity waves. The attempts to show how this 
can affect the composition changes in the lower atmosphere have not been 
successful. In view of these difficulties it is reasonable to make 
specific assumptions about the composition changes at the lower boundary 
and study their affects on the structural parameters of the neutral and 
ionized gas in the thermospheric region. The validity of these assumptions 
may then be tested by comparing the theoretical calculations with the 
observational data. This is the approach adopted in this paper. 

Discussion of the Results 

The Figures la and lb show the diurnal variations in I^F^ and the total 
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electron content n t (integrated from 120 to 1000 km) for the three condi- 
tions Q, D.| and D^. The corresponding variations in the neutral density 
at 200 km and the exospheric neutral temperature are shown in Figures 2a 
and 2b. The general inferences which can be drawn from these plots are as 
follows : 

1. With the boundary conditions specified in Table 1, both N F« 

m2 

and n t have decreased for the and conditions, compared to the Q 

condition. The decrease in (about 40%) is comparatively larger than 

the corresponding decrease in n^ (about 20%). A slight discontinuity in 
is caused by the development of the G condition to be discussed later. 

2. Both the neutral density and the exospheric temperature T^ 
have increased considerably during the condition. On the other hand 

the increase in the exospheric temperature during the condition is rather 

small and the density in fact is slightly lower than for the Q condition. 

It is important to note that T» and y at 200 km for the Q and 

conditions are substantially the same and are considerably different from 

the D 0 condition. Nonetheless the N F 0 and n. values for the D, and D« 

2 m 2 t 12 

conditions are about the same and are quite different from N m F 2 and n^ for 
the Q condition. 

These results merely point towards the difficulty of constructing 
ionospheric and atmospheric models based only on the neutral density as 
inferred from satellite drag measurements. If the density shown in Figure 2a 
would have been measured by satellite drag, the result would have been that 
the Q and D-j measurements refer to about the same atmospheric conditions, 
while the D-j and measurements refer to different atmospheric conditions. 
The conclusion then would have been (since all the other parameters are 
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assumed to be equal) that the ionospheres corresponding to Q and are 
the same and are different from the ionosphere corresponding to D^. As 
we see here, just the opposite is the case. This indicates that the total 
neutral density does not give enough information to be the basis of at- 
mospheric models. These points will be elaborated further in other examples. 

The ion composition in the altitude range of 120 - 400 km for the 
three cases (Q, D-| and D ? ) is shown in Figures 3-5. Each figure gives the 
altitude profile of 0 , N0 + along with the total electron density 

(sum of 0 + , NO and H + ) for 9.00, 12.00 and 15.00 hours local time 
on a linear scale. Figure 3 illustrates the development of the normal 
ionospheric layers F-j and F^ during the day. The electron density 
near the F 2 peak gradually increases towards the afternoon as a result 
of the increase in 0 + . The formation of the F-j layer is almost entirely 
determined by N0 + which attains its maximum values near noon. From 
Figures 4 and 5 which in our calculations represent stormtime conditions, 
we note that the 0 + layer is considerably depressed, causing a depression 
in the electron density near the F 2 peak. 

The development of the various ionic layers over the day during the 
D-j and D 2 conditions produces a rather interesting sequence of electron 
density profiles. With a single F layer in the morning hours to two layers 
(F-j and F 2 ) in the noon and afternoon hours, it shows the development of 
a pattern which has been commonly known as a G condition during magnetic 
storms. The G condition results from a situation in which the critical 
frequency of the F 2 "layer falls below that of F-j . A ground based ionosonde 
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in this case records only the reflection from the F^-layer. The situation 
is quite reversed when "seen" from the topside ionosphere. The ionogram 
in this case will show both F-j and F^ traces - a situation never encountered 
in the normal condition. A study of the G condition from topside ionograms 
has recently been made by Herzberg et al. (1969) and fatkullin (1970). The 
calculations presented in Figures 4 and 5 illustrate the physical situation 
for the development of the G condition and the recovery from it. 

The changes in ion composition during disturbed conditions in the 
altitude range below 300 km are not known experimentally. The results pre- 
sented here show that the most pronounced effect occurs in the 0 + layer 
which is deputed considerably during storm conditions. This is accom- 
panied by an increase in the N0 + concentration, resulting in a slight in- 
crease in electron density near the F^-peak. This is shown in Figure 6 
where the ion composition for the Q and conditions is presented for 15:00 
hours. The condition is not significantly different from and has 
been omitted in Figure 6 for clarity. 

Figure 7 shows the altitude profiles o* the neutral density and the 
ratio of [0] to [N^] for the three cases 0, D-j and 0^ corresponding to 15:00 

hours. The exospheric temperatures T for the three cases as obtained 

00 

from the solution of the heat conduction equation are also shown appropri- 
ately. It should be emphasized here again that all the parameters for the 
three cases are identical except for the differences in the assumed values 
of [0], [0 2 ] and [N 2 ] at the lower boundary. The D-j and Q conditions which 
differ only with respect to the atomic oxygen density at 120 km show no 


v 
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significant change in the neutral density up to 300 km. The difference 
in the exospheric temperature is only about 40°K. The ionospheric models 
constructed from the two conditions, however, are significantly different. 
The comparison of and presents the opposite situation. The two 
cases lead to completely different atmospheres. Their exospheric temper- 
atures are respectively 1084 and 1183°K. However, the ionospheric models 
constructed from these two cases are very similar. The understanding of 
this apparent dilemma is obtained by comparing the altitude profiles of 
[0]/ [N^] for the three cases. We note from Figure 7 that at any given 
altitude [0]/ [N^] is much larger for Q than for D-j or D^. 

The series of events which follow a change in [0]/[N 2 ] at the lower 
boundary is very complicated to describe. A certain insight, however, can 
be obtained from the following discussion: The main source of the thermal 

energy, responsible for heating the neutral atmosphere, is provided by the 
solar ultraviolet radiation. The U.V. is absorbed by the neutral consti- 
tuents of the lower atmosphere, consisting mainly of atomic and molecular 
oxygen arid molecular nitrogen. For a given U.V. flux the generation of 
thermal energy depends only on the relative concentration o* 0 2> and 0. 
Out of these three constituents atomic oxygen also acts as a heat sink 
losing energy via infrared radiation. Since the maximum thermal energy 
is generated in the altitude range around 120 km where N 2 is the dominant 
constituent, the increase in [N^] and (or) decrease in [0] results in a 
net increase of the thermal energy. This gives rise to a corresponding 
increase in the thermospheric temperature which, due to heat conduction. 
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extends to high altitudes. The three neutral constituents 0, and 0 2> 
through rapid diffusion, almost instantaneously attain a new diffusive 
equilibrium distribution, corresponding to this increase in temperature. 

The net result then is an increase in 0^ and throughout the altitude 
region under consideration. In the case of 0 this increase is reflected 
only in the high altitude region. A larger value of [N^] and [0^] and a 
smaller value of [0] results in less production and more loss for 0 + 
causing depletion of this layer. Since N0 + is produced at the expense of 
0 + , an increase in the loss rate for 0 + will cause an increase in N0 + . 

SUMMARY AND CONCLUSION 

Through the simultaneous solutions of a system of basic ionospheric 
and atmospheric equations we were able to simulate many of the observed 
features of the ionosphere for quiet and disturbed conditions. The main 
findings of this paper can be summarized as follows: 

1. A decrease in [0]/[N 2 ] at the lower boundary (120 km), arising 
from the decrease of the turbopause level during a magnetic storm, triggers 
a complex series of events resulting in the depletion of the 0 + layer 
and enhancement of N0 + . The development of these layers during the day can 
occasionally produce an interesting sequence of electron density profiles 
giving rise to the so called G condition. 

2. The ionospheric parameters are much more sensitive to changes in 
[0]/[N 2 ] than to changes in the absolute values of [0] and [N 2 ]. In the 
case of the neutral density, the situation is reversed. Since the iono- 
spheric parameters can be measured much more accurately than the atmospheric 
parameters, it is very useful to consider their properties i,i constructing 
the atmospheric models. 
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3. Because of the variability of the turbopause level during a 
magnetic storm, it is not possible to compute the neutral temperature 
from the density measurements alone. The temperatures derived on the 
assumption of the invariant boundary cannot be considered to be a measure 
of the gas temperature and as such have very little physical significance. 

The examples we chose in this paper to describe the stormtime behavior 
of the ionsophere correspond to what is generally called the negative 

t 

phase of the storm, a phenomenon more common in middle and high latitudes. 
In low latitudes, the storm is generally accompanied by an increase in 
f o F^, giving rise to the positive phase. If the stormtime changes in the 
ionosphere are caused by the movement of the turbopause level, it is easy 
to visualize that the positive and negative phase of the storm result from 
the rise and fall of this le vel caused by the global circulation. 

In this paper we limited our discussion to the changes in the neutral 
and ion compositions in the altitude range of 120 - 400 km. Even though 
the solutions of the heat conduction equations for electrons and ions (0 + , 
H + ) were an integral part of the solutions of the system of equations used 
in our study, we did not present any discussion on the behavior of electron 
and ion temperatures during magnetic storms. Our calculations showed that 
these temperatures increased during the 0^ and conditions following the 
pattern of the exospheric neutral temperatures. As we pointed out earlier 
in this paper, there are two principal mechanisms which are effective in 
depositing energy in the ionosphei^ during magnetic storms. These are: 

(1) The Joule dissipation of the ionospheric currents and; (2) the thermal 
conduction downward from the magnetosphere to the ionosphere '"'ong the geo- 
magnetic field lines. In this paper we investigated only the consequences 
of the first mechanism. The second mechanism is responsible for heating 
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the electron gas more efficiently, occasionally increasing its temperature 
to a point wher. it can give rise to a red arc condition through collisional 
excitation of atomic oxygen. Because of tne general interest in the problem 
of the mid-latitude arc, the mechanism for heating the electron-ion gas 
merits a separate consideration. This will be the topic of a second paper. 
We should like to conclude this paper, however, with the statement that 
the ion and neutral composition changes discussed here arise primarily 
from the changes f.n the turbopause level. The increase in the electron 
and ion temperatures as a result of the downward conduction have practi- 
cally no influence on the changes in ion and neutral composition helow 
about 400 km. 
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FIGURE CAPTIONS 


Figure la 

Figure lb 

Figure 2a 

Figure 2b 

Figure 3 
Figure 4 

Figure 5 

Figure 6 

Figure 7 


Diurnal variation of , the electron density at the F-peak, 
for the cases Q, D 1 and 

Diurnal variation of the columnar electron content up to 
1000 km for the cases Q, D-j and D^. 

Diurnal variation of the neutral density at 200 km for the 
cases Q, D^ and D^ 

Diurnal variation of the exospheric neutral temperature for 
the cases Q, D 1 and D 2< 

Changes in the ion composition during the day for the Q-case. 
Changes in the ion composition during the day for the D-j-case, 
showing the development of a G-condition around noon. 

Changes in the ion composition during the day for the D^-case, 
showing the development of a G-condition around noon. 
Comparison of the ion densities at 15 hours local time for 
the cases Q and D^. 

The neutral density and the ratio [0]/ [N^D as a function of 
height at 15 hours local time for the cases Q, D-j and D^. 
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